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Abstract. Aerodynamic ground effect in flapping-wing insect flight is of importance 
to comparative morphologies and of interest to the micro-air-vehicle (MAV) 
community. Recent studies, however, show apparently contradictory results of either 
some significant extra lift or power savings, or zero ground effect. Here we present a 
numerical study of fruitfly sized insect takeoff with a specific focus on the significance 
of leg thrust and wing kinematics. Flapping-wing takeoff is studied using numerical 
modelling and high performance computing. The aerodynamic forces are calculated 
using a three-dimensional Navier-Stokes solver based on a pseudo-spectral method with 
volume penalization. It is coupled with a flight dynamics solver that accounts for the 
body weight, inertia and the leg thrust, while only having two degrees of freedom: the 
vertical and the longitudinal horizontal displacement. The natural voluntary takeoff 
of a fruitfly is considered as reference. The parameters of the model are then varied to 
explore possible effects of interaction between the flapping-wing model and the ground 
plane. These modified takeoffs include cases with decreased leg thrust parameter, 
and/or with periodic wing kinematics, constant body pitch angle. The results show 
that the ground effect during natural voluntary takeoff is negligible. In the modified 
takeoffs, when the rate of climb is slow, the difference in the aerodynamic forces due to 
the interaction with the ground is up to 6%. Surprisingly, depending on the kinematics, 
the difference is either positive or negative, in contrast to the intuition based on 
the helicopter theory, which suggests positive excess lift. This effect is attributed 
to unsteady wing-wake interactions. A similar effect is found during hovering. 
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1. Introduction 

The aerodynamic forces of an air vehicle or an animal may be affected by the gronnd 
proximity. This phenomenon, known as the gronnd effect, has been extensively stndied 
for aircraft [3] and rotorcraft j20]. Althongh the effect varies depending on many design 
parameters, the general trend is an increase in lift and pitching moment, and a decrease 
in drag. The effect decays as the distance from the gronnd increases, and vanishes at 
a distance slightly larger than the characteristic length of the vehicle. For example, for 
a hovering helicopter, the excess thrnst vanishes if the distance to the gronnd exceeds 
1.25 times the diameter of the main rotor izg. 

Rayner [25] proposed a hxed wing lifting line theory for forward flight of birds, bats 
and insects. His analysis snggested that flight in gronnd effect provides performance 
improvements, if the flight speed is not too low. However, this theory conld not 
be applied to hovering or slow forward flight at very low height, since it neglected 
flapping motion. Normal hovering in gronnd effect was considered by Gao and Ln 
|13] . They carried ont two-dimensional nnmerical simnlations of hovering and identihed 
three regimes: force enhancement, force rednction, and force recovery, depending on the 
distance from the gronnd. Lin et al. considered clap-and-fling near the gronnd and 
fonnd force enhancement at all distances. A three-dimensional nnmerical simnlation 
of frnitfly hovering was carried ont by Maeda and Lin [23]. An increase in lift and a 
rednction in power was fonnd. A signihcant vertical force was generated on the insect’s 
body dne to the ‘fonntain effect’. Energetic savings have also been reported for a 
hnmmingbird hovering in gronnd effect HZI. 

Several stndies considered pitching-plnnging foils near a solid wall or a free snrface 
[28] l2T] l32] . This conhgnration is relevant to hsh swimming as well as forward flapping 
flight. The gronnd effect mainly consists in enhanced propnlsive force. However, it also 
generates a non-zero vertical force dne to asymmetry. 

The main motivation for this stndy comes from the fact that the gronnd proximity 
is natnral for takeoff and landing. These manoenvres, nnlike hovering or forward flight, 
are characterized by gradnal change of distance to the gronnd. The ‘dynamic’ gronnd 
effect in these circnmstances may be different from the ‘static’ effect at a constant 
distance US- This difference may be even larger for flapping wings than for hxed wings, 
becanse animals vary their wing kinematics dnring takeoff. 

So far, the gronnd effect dnring takeoff has been assessed for very few insects only. 
It was fonnd negligible for bntterhies {Pieris rapae [1], Papilio xuthus [22]), a dronehy 
{Eristarlis tenax) [7], and a frnithy {Drosophila virilis) [6], bnt signihcant for a beetle 
{Trypoxylus dichotomus) [30]. The disparity can be attribnted to signihcant diherences 
in the size, morphology and kinematics of these insects. Thns, onr work is motivated by 
the apparently contradictory conclnsions on the signihcance of the gronnd ehect that 
conld be fonnd in the animal hight literatnre. It is important to identify the parameters 
that make the gronnd ehect strong or negligible. 

In the present stndy, we consider a nnmerical model having the morphology of a 
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fruitfly, with variable wing kinematics and leg parameters. Our objective is to determine 
if the ground effect can be signihcant for this model, and which conditions can lead 
to it. We thus explore the parameter space of the model and perform numerous 
numerical simulations using FluSi [10], which is an open source software available 
on https://github.com/pseudospectators/FLUSI/tree/plos_one_ground_effect, 
First, for completeness, we revisit the voluntary takeoff of a fruitfly analyzed in [6]. 
The main difference with respect to [6] is the use of a flight dynamics solver. We then 
compare takeoffs with modihed parameters of the leg thrust model and wing kinematics. 
Finally, we consider hovering as a limiting case of very slow takeoff. 

The paper is organized as follows. In section |2l we describe our computational 
approach and the takeoff parameters used in this study. The results are presented in 
section m hrst for a natural voluntary takeoff, then for modihed takeoffs and for hovering 
hight. The main conclusions are summarized in section IH 

2. Methods 

2.1. Morphology and kinematics 

In this work, we consider a fruithy having mass m = 1.2 mg and wing length R = 2.83 
mm, which are the values reported by Chen and Sun [6]. The body is modelled as a rigid 
solid, and the wings are modelled as rigid hat plates. This approximation is accurate 
for drosophila during voluntary takeoh mu, though it occasionally fails during fast 
escape manoeuvres [12]. The wing contour used in this study is shown in Fig. [U^a). It is 
adapted from [6]. Its mean chord length is equal to c = 0.85 mm. The body is generated 
by sweeping a circular section of variable radius along a curvilinear centreline (an arc). 
The body has approximately the same dimensions as in [6] . The side view of the body 
is shown in Fig. HI) 6). Even though the yaw and roll angles can eventually become 
large during takeoh, there is no signihcant trend for all takeohs. Hence, to simplify the 
problem, we assume bilateral symmetry. Therefore, the body orientation is fully dehned 
by the pitch angle f between the body and the horizontal axis, see Fig. [I](c). The wing 
kinematics is described by three angles: 0, a and 9, measured with respect to the stroke 
plane, as shown in Fig. [I](d). The positional angle f dehnes the motion of the wing tip 
projection on the stroke plane. The deviation (elevation) angle 6 dehnes the deviation 
of the wing tip from the stroke plane. The feathering angle a dehnes the rotation about 
the longitudinal axis of the wing, and it is related to the geometrical angle of attack 
{AoA) as a = 90° — AoA during downstroke and as a = 90° + AoA during upstroke. It 
is convenient to refer to an ‘anatomical’ stroke plane angle r], i.e., to assume that the 
inclination of the stroke plane against the body axis is held at a constant angle for any 
motion of the body. 

Since the main focus of this study is the ground ehect, it is important to ensure 
that the time evolution of the distance to the ground is consistent with the forces acting 
on the insect. For this reason, in our computations, unlike in [6], the position of the 


Aerodynamic ground effect in fruitfly sized insect takeoff 


4 





Figure 1. Schematic drawing of the morphological model, (a) Wing contour. 
Coordinates are normalized to the wing length R. (b) The body is generated by 
circular sections of variable radius, which changes depending on the position along the 
centre line (dash-dotted arc). The body axis ObXb is the thorax-abdomen principal 
axis, approximately. The coordinates of the wing pivot points in the body frame of 
reference ObXbytZb are (—0.07i?, ±0.18i?, 0.115i?). Note that, since P{t) is prescribed 
in the simulations and only (xc, z^) are dynamically calculated, the position of the 
body point of reference with respect to the body contour is chosen arbitrarily, (c) The 
insect’s position with respect to the ground is described by the body point of reference 
coordinates {Xc,Zc) and the position angle (3. (d) Definition of the wing’s angles with 
respect to the stroke plane frame of reference OspXspVstZsp- The origin Osp is the wing 
pivot point. 


insect is dynamically computed as opposed to be prescribed. We compute the position 
of the body point of reference (xc, zf), see Fig. [I](6), from Newton’s 2nd Law, 


m- 


dfxc 


— Fax + Ft 


d2; 


ixy 


m- 


dF 


= Faz + Fi^ - mg, 


( 1 ) 


where {Fax, Faz) is the aerodynamic force, {Fix, Fiz) is the leg thrust, subscripts x 
and correspond, respectively, to the horizontal and vertical components, m is the 
insect’s mass and g is the gravitational acceleration. Equations o are integrated 
using the adaptive second order Adams-Bashforth scheme [26], simultaneously with 
the incompressible Navier-Stokes equations. We dehned the positive 2 ; direction to be 
upwards and the positive x direction to be forwards (see hgureHJc). 
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The aerodynamic forces Fax and Faz are obtained by solving the three-dimensional 
incompressible Navier-Stokes equations. The no-slip boundary condition at the body 
and wings surfaces is imposed using the volume penalization method [1], and the 
penalized equations are solved using a classical Fourier pseudo-spectral method. More 
details about the solver and the generic insect model, including a numerical validation 
case of fruitfly hovering, can be found in [10]. Numerical validation of the ground plane 


modelling using the volume penalization method is described in Appendix A 


Ly and L 


The computational domain in the present study is a rectangular box with sides Lj., 
Suitable values of Lx, Ly and L^, in terms of accuracy and computational 
efficiency, depend on the motion of the insect within the domain. Therefore, different 
values are used in different simulations, as described later in the text. The domain is 
discretized using a uniform Cartesian grid. Periodic boundary conditions are applied 
on all sides of the domain, as required by the Fourier discretization. Vorticity sponge 
boundary conditions are imposed at the left, right, rear and front sides of the domain, 
as explained in mm, in order to minimize the effect of the finite domain size. The 
ground surface is modelled as a solid layer at the bottom of the domain, which by 
periodicity also imposes the no-slip on the top of the fluid domain. We have carried 
out numerical experiments to ensure that, in the numerical simulations presented in 
this paper, the domain size is sufficiently large, i.e., its further increase does not change 
significally the forces. The dimensions that we chose are also comparable with the size 
of the mineral oil tank used in the experiments with a mechanical model [9]. 


2.3. Leg thrust 

The model of the leg thrust employed in the present study is a slight modification of 
the compression spring model proposed in [1]. We assume that takeoff begins from rest 
and starts at time t = tg, which can be estimated from the initiation of the legs motion 
in the video sequences shown in [6j. The two components of the force are given by 

Fix = Z cot fi, Fi^ = Z. (2) 

The magnitude of the leg force is assumed to depend on the vertical component of the 
leg extension ( = zdt) —Zcitf) only. The force is supposed to be distributed between the 
three pairs of legs such that its change with horizontal displacement can be neglected. 


Z = 


{L,-C)K, forC<T,, 

0 for (■ > L^, 


(3) 


where Li is the maximum leg extension length, i.e., the difference between the values 
of when the legs are fully extended at takeoff and when the insect is at rest. When 
the legs are fully extended, C, = Li, the legs lose contact with the ground and the force 
drops to zero. This length is estimated using video sequences in [6] to be equal to 
L^ = 1.24 mm. The spring stiffness varies in time: it increases from Kf before 
takeoff to iF/ after takeoff. The initial value Kf = mglLi ensures that the insect is 
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in equilibrium before takeoff, when the aerodynamic force is zero. The final value 
is a parameter of the model that controls the maximum leg thrust. Its value can be 
estimated from the climb velocity at the beginning of takeoff, shown in, e.g., [6]. It may 
also be estimated from jumps of wingless flies [33l |5] for a slightly different fruitfly, D. 
melanogaster. We assume the time evolution of of the form 

{ Kf for t <ti, 

K+ — K~ 

Kf - - - —{t-tf) ioi ti <t <t(, +n, (4) 

Te 

for t > t£ + Ti- 

The transition time ti can be equal to zero, in which case the leg force increases 
impulsively at the beginning of takeoff. However, measurements of the leg force [33] 
suggest a gradual increase which can be accounted for by setting ti to a value larger 
than zero. The value ti = 1.3ms results in the gradient dF^/dt consistent with the 
experimental data shown in [33|. The direction cfi also changes in time. Before takeoff, 
when the insect is at rest, the force is applied only in the vertical direction, i.e., ff = 90°. 
During takeoff, the horizontal component is non-zero, in general. We assume a time 
evolution of the form 

{ ff for t < ti, 

4>f F— - -{t-tfj ioT te <t <te + Te, (5) 

Te 

ff for t >te + Te- 

The values of the leg thrust model parameters used in our numerical simulations are 
given in table [H 

3. Results and Discussion 

3.1. Summary of the numerical simulations 

The starting point for our study is the voluntary takeoff, as it is shown in section 13.21 
(in agreement with [6]) that the ground effect is very small in that case. It is much 
smaller than during hovering (cf. |23|). We conjecture that this difference is due to 
the large takeoff vertical velocity, which is mainly the result of the leg thrust. To 
test this hypothesis, in section 13.31 we discuss a situation in which the legs produce 
less force and the insect takes off slower. The ground effect becomes signihcant. The 
vertical force increases during the first two wingbeats due to the ground effect, but 
slightly decreases later on. We then carry out a parametric study using periodic wing 
kinematics in section 13.41 and hnd an even stronger adverse ground effect. Finally, 
in section 13.51 we find similar trends during the first wingbeats in hovering flight, 
which can is considered as a limiting case of takeoff with zero rate of climb. Table [T] 
summarizes the parameters of the different cases considered in the present study. 
Datasets for the ‘voluntary’ and ‘simplified’ cases can be downloaded from ResearchGate 
at http://dx.doi.org//10.13140/RG.2.1.2145.2562, 
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Table 1. Parameters of the takeoffs considered in the present study. 


Case 

Kinematics 

Leg model 

Name 

(j), a, 9 

13 

vf 

Zc{0), mm 

Li, mm 

K+, N/m 


ti, ms 

Ti, ms 

Voluntary 

figja 

figl2] 

62 

1.08 

1.24 

0.165 

84 

4.2 

1.3 

Slow 

figja 

figl2] 

62 

1.08 

1.24 

0.041 

84 

4.2 

1.3 

Simplified 

figlZKa) 

46.3° 

32 

3.11 

1.24 

0.0095...0.043 

84 

0 

1.3 

Hovering 

^gWLa,h) 

55° 

55 

2.07 

N/A 

N/A 

N/A 

N/A 

N/A 

Case 

Numerical parameters 


Name 

Lx 

Ly 

L. 

N. 

Ny 

N, 

e 

Voluntary 

5R 

5R 

8R 

640 

640 

1280 

2.5 • 10-^ 

Slow 

5R 

5R 

6R 

640 

640 

768 

2.5 • 10-^ 

Simplified 

4R 

4R 

6R 

512 

512 

768 

2.5 • 10-^ 

Hovering 

8R 

8R 

4R 

864 

864 

432 

2.5 ■ 10“^ 


3.2. Voluntary takeoff 


In this section, we consider voluntary takeoff of a fruitfly with the parameters as in the 
first line in Table [H This case shows some important general features of fruitfly takeoff 
such as the first wingbeat cycles beginning while the legs extend. Therefore it is likely 
that, despite some variability in voluntary takeoffs, the ground effect in general remains 
of the same order of magnitude in natural circumstances. 

The values of the body and wing angles are taken from one of the cases documented 
in [6]. However, the wing motion in [6] is not exactly symmetric. Therefore, the time 
series of 0, a and 6 that we use for both wings correspond to the left wing data shown in 
[6]. Fig. [2] presents the time evolution of the wing positional angle 4>{t), the feathering 
angle a(t), the elevation angle Oft) and the body pitch angle (3ft), which are prescribed 
in our numerical simulations. The angle between the horizontal plane and the stroke 
plane, r] — (3, is also shown for reference. 

Even though the wing motion is not exactly periodic, it is useful to introduce the 
wing beat frequency. When calculated using the average wing beat cycle period over 
the hve cycles shown in Fig. [2l it is equal to / = 169 Hz. Similarly, the average 
wing beat amplitude is equal to $ = 134°, and the characteristic wing tip velocity is 
U = 2^Rf = 2.23 m/s. The kinematic viscosity of air, equal to z/ = 1.45 ■ 10“^ m^/s 
yields the Reynolds number Re = Uclu = 131. Note that U and Re do not account for 
the forward speed of the body. 

The computational domain size is equal to = Ly = 5R, = 8R, where R is the 

wing length. The influence of the domain size in the vertical direction is discussed in 


Appendix B[ The number of grid points in each direction, respectively, is = Ny = 640 
and Nz = 1280. The penalization parameter is e = 2.5 • 10“^ (for details see, e.g., |lljb 
The aerodynamic ground effect is evaluated by comparing two numerical 
simulations with two different values of the initial distance from the body point of 
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Figure 2. Time evolution of the angular position of the body and of 
the wings during the voluntary takeoff. Gray shaded regions correspond to 
downstrokes. 77 —/3 is the angle between the horizontal plane and the stroke plane, i.e., 
the global stroke plane angle [ 2 ^ . 


reference to the ground: Zc{0) = 0.38i? and 2R, which we denote ‘in ground effect’ 
(IGE) and ‘out of ground effect’ (OGE), respectively. The hrst case corresponds to 
a takeoff from a flat ground surface, with Zc{0) being consistent with the data in [6]. 
In the second case, the leg model behaves as during takeoff from the ground, but the 
aerodynamic interaction between the insect and the ground is weak because of the large 
distance. This case may be interpreted as takeoff from a perch that provides enough 
support for the legs but has a small surface, such that the aerodynamic interactions 
are negligible. With the distance equal to 2R or greater, the ground effect is negligible 
during hovering [22]. The circulation of the wake vortices is mainly determined by the 
integral aerodynamic force, therefore it is not larger during takeoff than during hovering, 
and the spatial rate of decay of the induced velocity is the same. Hence, the ground 
effect with the distance equal to 2R is likely to be negligible during takeoff. The influence 
of the ground on the shape of the vortices is only visible during the 2nd wingbeat and 
later on. This difference is localized to the vicinity of the ground plane. Since the insect 
is relatively far from the ground by that time, this difference is unlikely to have any 
influence on the aerodynamic forces. 

Fig. ma) shows the fruitfly model and the wake, IGE and OGE, at 4 subsequent 
time instants. The vortices created by the wings and the body are identihed as the 
volume of fluid enclosed by the iso-surfaces of the Q-criterion. At f = 0, the air is at 
rest. The insect body is almost horizontal. The wings are in a pre-takeoff position from 
which they begin the hrst downstroke after t = 4.1 ms. The time t = 9.2 ms corresponds 
to the hrst reversal from downstroke to upstroke. Because of the small body pitch angle 
/3, the stroke plane is ehectively vertical. In addition, the wing tip speed during the 
hrst downstroke is smaller than during all subsequent strokes. Therefore, the vertical 
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(a) 
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plane 


t = 9.2 ms, 1st wingbeat cycle 




ground Dlanp 


t = 12.8 ms, 2nd wingbeat cycle 
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ground Dianp 
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(e) 
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Figure 3. Voluntary takeoff, (a) Visualization of the wings, body and ground 
surface, and the wake at 4 subsequent time instants. Blue semi-transparent iso-surfaces 
show the Q-criterion, QjP = 15. (6) Vertical and horizontal displacement. To obtain 
distance Zc from the ground for the IGE case, add 1.08 mm. The black dash-dotted 
line indicates = R. (c) Components of the leg force, (d) horizontal and (e) vertical 
components of the aerodynamic force and (/) the aerodynamic power. The black dash- 
dotted line in Fig. (e) indicates the weight. Solid circles connected by dotted lines 
show wingbeat cycle averages. The results for OGE and IGE are shown, but the curves 
in Fig. (6-/) overlap because the difference is negligible. 
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aerodynamic force is small, but the body lifts noticeably because of the leg thrust. The 
time t = 12.8 ms corresponds to the second upstroke. At this point, the distance from 
the body point of reference to the ground Zc is already larger than the wing length R. 
Therefore, the aerodynamic interference with the ground is expected to be very small. 
Note that the kinematics during the first two wingbeat cycles are a transient. After 
that, the time evolution of the wing angles approaches a periodic regime and the stroke 
plane becomes less inclined with respect to the ground, see Fig. [2l 

The displacement of the body point of reference is shown in Fig. [3](6). It presents 
the evolution of the vertical component f{t) = Zc{t)—Zc{0) and the horizontal component 
f{t) = Xc{t) — Xc{0) over time for the cases IGF and OGE. The curves overlap. In both 
cases, at the end of the 4th wingbeat cycle, t = 28.4 ms, the insect gains 8.7 mm of 
altitude and propels 1.6 mm forward. These numbers are consistent with the trajectories 
shown in [6]. The displacements IGE and OGE differ by less than 1%. Therefore, the 
ground effect on ( and f is indeed negligible. 

Fig. Et^c) shows the two components of the leg force. At f = 0, the vertical 
component of the leg force is equal to the weight and the horizontal component is 
zero. The jump is triggered at ti = 4.2 ms. At time te + T£ = 5.5 ms, both components 
reach their peaks. After that the force decreases and vanishes at t = 9.3 ms, when the 
legs lose contact with the ground. Note that the leg thrust can, in principle, be different 
for the takeoffs IGE and OGE, because the leg model depends on the aerodynamic force 
via Zc{t). However, for the voluntary takeoff considered here, there is no influence of 
the ground effect. 

The vertical and the horizontal components of the aerodynamic force are shown 
in Fig. [3](d) and (e), respectively. Over the first four wingbeat cycles, the wingbeat 
averaged aerodynamic forces are significantly lower than the weight. This can be 
explained by the large initial rate of climb due to the leg thrust, which cannot be 
supported by the wings. Even during the fourth wingbeat, the wing force is equal to 
29% of the weight. The vertical acceleration is therefore negative after the legs lose 
contact with the ground, and the rate of climb slowly decreases. The ground effect is, 
again, negligible. Even during the first wingbeat cycle, when the wings approach the 
ground surface, the difference in the instantaneous vertical force between IGE and OGE 
is at most 0.0005 mN, i.e., about 4% of the weight. The wingbeat cycle averaged forces 
differ by less than 1% of the weight. 

Fig. [3](/) displays the time evolution of the aerodynamic power, when operating 
IGE and OGE. Note that, in this study, we do not consider the inertial power because 
the wings have the same kinematics in both cases, IGE and OGE. Therefore, the inertial 
power is the same. The aerodynamic power is the aerodynamic component of the power 
required to actuate the wings, 

P = -Mi-{ni-ftb)-Mr-{nr-nb). ( 6 ) 

In ([6|), Ml and Mr are the aerodynamic moments of the left and of the right wing, 
respectively, relative to the corresponding pivot point. O; and are the angular 
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velocities of the wings and Qb is the angular velocity of the body. All vectors are 
taken in the laboratory frame of reference. P is positive if power is consumed. We 
hnd that it is positive during most part of the takeoff (see Fig. |3]f). Only at the 
reversals during the first two cycles, when the body velocity is still small, P is slightly 
negative. During the 2nd wingbeat, the mean body-mass specihc aerodynamic power 
is equal to Pf = Pave/^n = 21 W/kg. Assuming that the muscles contribute to 
30% of the body mass, the mean muscle-mass specihc aerodynamic power is equal 
to Pf = Pave/{0.3m) = 69 W/kg. The relative difference in the cycle averaged values 
between IGE and OGE is less than 0.5%. 

We conclude that the ground effect is unimportant for the voluntary takeoff, a 
result which is in agreement with [6] . This is mainly a consequence of rapid acceleration 
during the hrst wingbeat cycle, when the legs produce a large vertical force. The main 
question of the next section is whether this scenario changes if the takeoff is slower and 
the insect remains near the ground for a longer time. The rate of climb at the beginning 
of takeoff is controlled by the leg model stiffness coefficient K/', and the horizontal 
velocity is controlled by the leg angle ff. 

3.3. Slow takeoff 

This section describes a modihed takeoff with the leg thrust coefficient decreased to 
K/f = 0.041 N/m (see the second line in Tabled]). Smaller Kf results in less leg thrust 
and slower climb, compared to the natural voluntary takeoff. Therefore we refer to 
this case as a ‘slow takeoff’. Any further decrease of Kf would result in a very close 
approach of the insect to the ground surface and, ultimately, to a collision which would 
require special treatment. Modelling of such collisions could be an interesting topic for 
future research (see [3T] for a review of structural modelling of insect wings, including 
impact modelling). For a fruitfly, collisions between the wings and the ground may be 
undesirable because of the large wingbeat frequency and light wing structure. Therefore, 
K/' = 0.041 N/m is an interesting limiting case. 

In the present numerical simulations, the computational domain size in x and y 
directions, the discretization grid step size and the penalization parameter are the same 
as in the previous section. The domain size in the vertical direction 2 ; is reduced to OR 
because the insect gains much less altitude by the end of the simulation. 

Fig. Wfo) shows the displacement of the body point of reference. The rate of climb 
is about one third of its original value and the insect only gains 3.9 mm by the end 
of the 4th wingbeat cycle. This is just slightly larger than the wing length R (2.83 
mm). The displacement is slightly larger for IGE than for OGE in both directions, 
horizontal and vertical. The time evolution of leg thrust is given in Fig. Wib). There is 
no visible difference between the two cases. The peak of the vertical force is equal to 
0.051 mN, which is about four times less than in the original voluntary takeoff discussed 
in section 13.21 

The time evolution of the instantaneous aerodynamic force, shown in Fig. |4l(c,d), is 
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(a) (b) 




time, ms 


time, ms 



0 5 10 15 20 25 30 



0 5 10 15 20 25 30 


Figure 4. Slow takeoff, (a) Vertical and horizontal displacement. Dash-dotted line 
indicates Zc = i?- (b) Components of the leg force, (c) horizontal and (d) vertical 
components of the aerodynamic force. The black dash-dotted line in Fig. (d) indicates 
the weight. Solid circles connected by dotted lines show wingbeat cycle averages. Note 
that, in Fig. (b) and (c), the lines for ICE and OGE almost coincide. 


qualitatively similar to the voluntary takeoff case considered previously. The difference 
between the cases OGE and IGE is negligible for the horizontal force (Fig. |4]c), but 
for the vertical force it reaches values as large as 0.0027 mN, i.e., 23% of the weight 
(Fig. md). Fig. [5] shows the difference between the wingbeat averaged forces in the 
cases IGE and OGE, normalized by the weight. The vertical force difference is shown 
in Fig. [5](fe). During the 1st wingbeat, the ground effect makes the total vertical force 
increase by almost 6% of the weight (red line). However, during the 2nd wingbeat, the 
extra force decreases to only 2% of the weight. During the 3rd, the 4th and the 5th 
wingbeats, the difference between the vertical forces IGE and OGE is very small and 
negative. The increase of the vertical force during the hrst wingbeat is mainly due to 
the wings (bine line). The extra force acting on the body is only abont 1% of the weight. 
However, later on, the contribntion of the body becomes important becanse it remains 
positive, whereas for the wings it becomes negative. The horizontal force difference, 
shown in Fig. [5](a), is positive, i.e., the propulsive force increases dne to the ground 
effect by about 2% of the weight, for all wingbeats. The contribntion of the body is np 
to 1% of the weight. For reference. Fig. [5] also shows the force differences dnring the 
volnntary takeoff. They are all smaller than 1%. 

The aerodynamic power, in the cases IGE and OGE, is compared in Fig. [6l The 
maximum difference is of abont 3% in magnitude for the slow takeoff, but less than 1% 
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0 5 10 15 20 25 30 0 5 10 15 20 25 30 

time, ms time, ms 

Figure 5. Slow and voluntary takeoffs. The difference between the cases IGE 
and OGE, in terms of the wingbeat cycle averaged aerodynamic force normalized by 
the body weight (a) horizontal and (&) vertical components. The forces acting on the 
wings and the body are shown separately. The total force, which is their sum, is also 
shown. 



0 5 10 15 20 25 30 

time, ms 

Figure 6. Slow and voluntary takeoffs. Aerodynamic power ratio IGE/OGE. 


for the voluntary takeoff. Considering the slow takeoff, the insect consumes more power 
when operating in ground effect (IGE) during the hrst two wingbeat cycles, but less 
power during the subsequent cycles. Overall, we hnd that the differences in the power 
are small. 

S.4- Takeoffs with simplified kinematics 

In the previous sections we noticed that the ground effect depends on the takeoff 
kinematics. We are mainly interested in the effects that might be generally applicable to 
a fruitfly sized insect takeoff. Therefore, in this section we consider parametric studies. 
They are performed using simplified periodic wing kinematics. The time evolution of 
the three wing angles over one wingbeat period, obtained by periodization of the last 
wingbeat in [6], is shown in Fig. Wfa). The wingbeat frequency is equal to / = 210 Hz. 
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(a) 



wingbeat time fraction 


- positional, (() 

- feathering, a 

- elevation, 6 

- 100% amplitude 

— — — 50% amplitude 


(b) 



wingbeat time fraction 


Figure 7. Wing kinematics, (a) Data obtained by periodization of the last 
wingbeat shown in [^. It is used in section 15^ Also, in section 1531 this kinematics is 
referred to as ‘PI’ (‘P’ for ‘Periodic’). The cycle begins from the upstroke. (6) Data 
adapted from [53]. In section 1331 it is referred to as ‘P2’. The cycle begins from the 
downstroke. 


This takeoff mode can be relevant to MAVs, for which the wing kinematics and the 
body angle during takeoff do not change as much as for the fruitfly (see, e.g, 0). The 
leg strength parameter is varied, resulting in a variation of the takeoff rate of climb 
Vt,o.- The body angle is constant and equal to (3 = 46.3°, the anatomical stroke plane 
angle is equal to rj = 32°. In these computations, we use = Ly = 4/2, = QR, 

= Ny = 512 and = 768, corresponding to more than 200 million grid points. The 
penalization parameter is equal to £ = 2.5 • 10“"^. 

Smaller Kf implies smaller rate of climb (Fig. |8]c) which leads to a more signihcant 
ground effect (Fig. (Hid). A striking feature of Fig. is a signihcant decrease of the 
vertical force during the 4th, 5th and 6th wingbeats, by up to 6%. The horizontal force 
varies slightly, by about 1% (see Fig. [8]a). 

The decrease in vertical force on the wings found in the IGF cases compared to 
OGE (Fig. IHf)) is an adverse ground effect which may be the result of complex wing- 
wake-ground interactions that depend on wing and body kinematics. Adverse ground 
effects have been reported previously for hxed-wing aircraft [16], but they consist in 
increased drag together with increased lift. However, for a two-dimensional ellipse with 
normal hovering kinematics PI, the mean vertical force decreases when the height from 
the ground is between approximately 1.5/7 and 4/7, where D is the chord length of the 
ellipse. For flapping wings, an adverse ground effect was found by Quinn et ah m- 
They considered an airfoil undergoing pitch oscillations in a closed-loop water channel 
with prescribed free-stream velocity. Such a conhguration represents a section of a bird 
wing in forward flight or a hsh £n. Even if the pitching motion was symmetric, the 
proximity of the ground broke the symmetry of the flow. Thus, the airfoil produced 
non-zero lift. The lift was positive if the distance to the ground was less than 40% of 
the chord length, but it became negative at larger distances, such that the lift force 
pulled the airfoil towards the ground. Nevertheless, the extra propulsive force due to 









Aerodynamic ground effect in fruitfly sized insect takeoff 


15 



Figure 8. TakeofTs with simplified kinematics, (a) Horizontal and (b) vertical 
difference between the wingbeat-averaged force IGE and OGE, normalized to the insect 
weight, (c) Vertical velocity of the body point of reference (rate of climb) versus time. 
(d) Maximum normalized force difference versus takeoff rate of climb at the moment 
when the legs lose contact with the ground. 


the ground effect was positive in all cases. Note, however, that the flows considered in 
[13] and [21] are effectively two-dimensional. 

In the present work, importantly, we hnd an adverse ground effect in a three- 
dimensional conhguration, which has not been previously recognized. The wingbeat 
cycle averaged vertical force of the wings in the ground effect is slightly larger during 
the first two cycles, but, as the insect flies away from the ground, the vertical force in 
the case IGE becomes less than that in the case OGE. 

Fig. [9] shows a comparison of the wake at the end of the 5th wingbeat cycle 
(t = 23.78 ms) in the 4 different takeoffs IGE with different values of Kf. It corresponds 
to the maximum decrease of the vertical force. There are noticeable differences between 
the vortices when the takeoff is fast and when it is slow. The part of the wake that 
approaches the ground deforms when it impinges on the ground. It then rolls up in a 
pair of vortex rings. Similar “ground vortices” are known in the context of helicopter 
rotor aerodynamics. In each of the 4 cases shown in Fig. [9l they have different strength 
and position with respect to the wings. Therefore, they induce the downwash of different 
strength. 

A detailed view of the flow near the wings is presented in figure [JO] It shows 
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{a) A:,+=0.0095N/m 


{b) A:i+=0.0184N/m 




(c) ^,^=0.0307N/m 


Ground 

vortex 


{d) ^r=0.0430N/m 


Figure 9. Flow visualization at the end of the 5th wingbeat (t = 23.78 ms) for 
the simplified kinematics cases. Iso-surfaces Q//^ = 15 are shown for 4 different 
takeoffs IGE: (a) Kf = 0.0095 N/m; (6) Kf = 0.0184 N/m; (c) Kf = 0.0307 N/m; 
{d) Kf = 0.0430 N/m. 


the pressure and the vorticity magnitude during the 6th wingbeat for two IGE cases 
with different Kf. Five time instants are visualized in Eve rows, respectively. The left 
column shows the pressure distribution over the surface of the insect, as well as over 
a semisphere of radius 0.9i/ centred the body reference point, for Kf = 0.0095 N/m. 
Figure fTOf ffl) is at f = 24.73 ms, during upstroke. The dark blue area near the leading 
edge that expands towards the wing tip is the trace of the leading-edge vortex (LEV), 
similar to the one discussed in [2]. The LEV at the downstroke is evident in hgure fTOf d) 
at f = 27.58 ms. The pressure distributions during the reversals are more complex 
fhgures flOh. c and e). 

The middle column of figure [lO] shows the pressure iso-contours sampled on the 
surface of the semisphere. Two cases are compared, Kf = 0.0095 N/m and 0.0430 N/m. 
This choice of Kf corresponds to the largest difference in the vertical force (see 
figure |8]&). For each of the cases, two isolines are drawn, p/pR^f"^ = —2 and —5. There 
is a large difference between the contours for different Kf in the far wake (wake far 
from the wings). However, in the near wake of a wing including the LEV, the difference 
is much smaller, which is consistent with the force changing by only a few per cent. 
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The iso-contours of the vorticity magnitude are compared in the right column. Here 
again, despite signihcant differences in the far wake, the LEV contours virtually overlap 
for K'^ = 0.0095 N/m and 0.0430 N/m. This shows that the ground effect has almost 
no influence on the dynamics of the LEV. 

3.5. Ground effect in hovering flight 

In this section, we simplify the kinematics even further. We consider hovering with 
the insect body being hxed. The flight dynamics solver is not used in this case. The 
distance from the body centre to the ground is equal to 0.48i? (where R is the wing 
length) for hovering in ground effect (IGE) and 2AR for hovering out of ground effect 
(OGE). The body pitch angle and the anatomical stroke plane angle are both constant 
and equal to 55°, such that the stroke plane is horizontal. The wing kinematics is the 
same as in the previous section, see Fig. Ha). We denote it as ‘PI’ kinematics. The 
wingbeat frequency is equal to / = 218 Hz. The hrst wingbeat starts from the upstroke, 
as done in [6]. 

In these numerical simulations we are interested in the long-time evolution of the 
aerodynamic forces, which after the initial transient eventually reach a periodic state. 
Most of the results known from the helicopter rotor theory [20] are obtained in reference 
to the periodic state, while the takeoffs considered in the previous sections of this paper 
(the slow takeoffs, in particular) last only for a few wingbeats. Therefore, it is instructive 
to consider the time evolution of the aerodynamic forces during hovering from f = 0 
until the time when the periodic state is reached. 

Since the time span of the numerical simulations presented in this section is 
large, it is necessary to increase the domain size in the horizontal directions. We set 
Lj. X Ly X Lz = 8R X 8R x 4/2, where z is the vertical direction. The number of grid 
points is Nx x Ny x = 864 x 864 x 432. The penalization parameter is equal to 
8 = 2.5- lO-P 

The quantity of interest is the ratio of the wingbeat averaged forces, IGE to OGE: 
Fz,ave,iGE/Fz^ave,OGE- This quantity is shown in Fig. [TTTai. The red solid line with 
“-I-” symbols corresponds to the total vertical force ratio. As already noticed in [23] . 
the vertical force during hovering in ground effect, Fzjge, reaches its periodic state 
signihcantly later than during hovering out of ground effect, Fz^oge- Therefore, the 
ratio of their wingbeat averages, Fz^ave,iGE/Fz^ave,OGE converges slowly with the number 
of wingbeats. It oscillates between 102% and 107%. After 27 wingbeats it reaches 106%. 

A similar comparison for the force generated by the wings is shown with a red 
dot-dashed line in Fig. fTTTffl). It was calculated by integration of the distributed forces 
over the wings only, in the same numerical simulations. Therefore, the aerodynamic 
interaction between the body and the wings is included. This force ratio drops from 
100.6% to 92.5% during the hrst 6 wingbeats, oscillates and then increases to 96.3%. 

The time evolution of the wake vortices generated by the insect is shown in Fig. [121 
The visualized time instants correspond to the end of the 1st, 2nd, ..., 25th downstroke. 
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Figure 10. Visualization of the leading-edge vortex during the 6th 
wingbeat for the simplified kinematics cases, IGE. Left column (a-e) shows the 
dimensionless pressure distribution over the surface of the model and over a semisphere 
of radius 0.9i? around the body point of reference. Middle column {f-j) shows the 
pressure iso-contours for two different takeoffs: Kf = 0.0095 N/m and 0.0430 N/m. 
Right column (fc-o) compares iso-contours of the dimensionless vorticity magnitude for 
the same two takeoffs. Time instants are t = 24.73, 25.68, 26.63, 27.58 and 28.53 ms 
{tf = 5.2, 5.4 5.6, 5.8 and 6, where / is the wingbeat frequency). 
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(a) 



Wingbeat cycles 


-1- PI, Total, 100% amp, 

— — PI, Wings, 100% amp. 

-1- PI, Total, 50% amp. 

— — PI, Wings, 50% amp. 

-1- P2, Total, 100% amp. 

— — P2, Wings, 100% amp. 

P2, Total, 50% amp. 
P2, Wings, 50% amp. 


(b) 



Figure 11. Hovering flight, (a) The ratio of wingbeat averaged vertical force IGE 
to OGE, Fz ave,iGE/Fz,ave,OGE- (^) The ratio of wingbeat averaged power IGE to 
OGE, PaveJGE/Pave,OGE- 


There are significant differences between the positions of the vortices over the first fonr 
time instants. The first wingbeat generates very strong vortex rings, that collide with 
the gronnd. Then they rebonnd dnring the second wingbeat, and parts of them moving 
npwards are still visible dnring the third wingbeat. The downwash prodnced by these 
vortices infinences the nearer wake dynamics and it is likely to be responsible for the 
decrease of the vertical force dnring the first few wingbeats. After the 10th wingbeat, 
the wake approaches its qnasi-periodic state. There are almost no visible differences 
between the visualizations at the end of the 20th wingbeat and at the end of the 25th 
wingbeat. 

The pair of numerical simulations (cases IGE and OGE) that we have discussed in 
the above paragraphs leads to the following conclusions. 

(i) Over the first 27 wingbeats, F^^aveyGE/Fz,ave,OGE varies within about 5% for the 
total force and 9% for the wings force. 

(ii) The wing generates less vertical force in the case IGE than in the case OGE (adverse 
ground effect). 

(iii) The body makes an important contribution to the total vertical force when 
operating IGE, which results in the excess total vertical force (positive ground 
effect). 

These conclusions are, of course, only valid for the particular wing shape and kinematics 
used in the simulation. Periodic flapping is only an approximation to the real insect 
wing motion, which varies from one wingbeat to another, and depends on many different 
conditions. To determine the effect of all existing fruitfly wing kinematics is beyond the 
reach of our numerical simulations. However, it is useful to compare a few different 


cases. 
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1st wingbeat 



1st wingbeat vortices 
approach the ground 




2nd wingbeat 1st wingbeat vortices 

rebound from the ground 



3rd wingbeat 


1st wingbeat vortices move 
far away from the ground 



4th wingbeat 





Figure 12. Time evolution of the wake during hovering. Iso-surfaces of the 
Q-criterion, Q/= 15, are shown at the end of the downstroke. ‘PI’ kinematics with 
100% wingbeat amplitude. 


We carried out numerical simulations with the wing kinematics used in [23] 
(abbreviated as ‘P2’ in the hgures). Note that, in this case, the first wingbeat begins 
from the downstroke, as shown in Fig. [ 71 (&). The wingbeat frequency is the same, 
/ = 218 Hz. The results of these numerical simulations are shown in Fig. [11] with 
green lines. They are qualitatively similar to the previously shown ‘PF case, but the 
values are systematically larger. F^^aveyGE/Fz^ave,OGE reaches 111% for the total force 
and 100% for the wings, such that there is no adverse ground effect after the periodic 
state is established. 

The adverse ground effect is rarely encountered in the aircraft or rotorcraft 
aerodynamics literature. However, in the context of flapping wings, it is not unusual. 
In the two-dimensional numerical simulations [13], a U-shape profile of the force ratio 
Fz,ave,iGE/Fz^ave,OGE versus k/c was found, where h is the distance from the wing centre 
to the ground and c is the wing chord. Fz,ave,iGE/Fz,ave,OGE was greater than 100% for 
hjc < 1.5, but less than 100% for hfc > 1.5, and the minimum ratio was of about 54%. 

In the three-dimensional model considered in the present paper, it is possible to 
partially reduce the three-dimensional effects by decreasing the wing beat amplitude. 
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One can then expect the adverse ground effect to be amplified. Indeed, this is what we 
hnd by decreasing the wingbeat amplitude by a factor of 2 (by rescaling the positional 
angle shown in Fig. [7] such that the amplitude of the positional angle is halved but the 
mean positional angle is unchanged, and the wingbeat frequency remains unchanged). 
The results are shown in Fig. [Hlwith blue lines. Now we have Fz^ave,iGE/Fz,ave,OGE < 
100% at intermediate times, for wings force and for the total force. The final periodic 
state produces a very slight excess of the total vertical force (less than 1%). Note that, 
in this reduced-amplitude case, the total force ratio decreases more than the wing force 
ratio. This indicates, not surprisingly, that the fountain effect becomes weaker when the 
wingbeat amplitude is reduced. Similar computations with the ‘P2’ kinematics show 
the same trend with an even larger decrease of F^^aveyGE/Fz^ave,OGE- 

The wingbeat averaged aerodynamic power ratio PavejGE/Pave,OGE is shown in 
Fig. [TIT M. Its variation is smaller than the variation of the force, and the computations 
suggest that its long-time limit is between 97% and 99%, in all cases that we have 
considered. The shape of the time evolution profiles of the power ratio is approximately 
similar to the time profiles of the wings vertical force ratio. This means that a local 
decrease of the vertical force ratio is accompanied by a decrease of the power ratio. 
Therefore, if the kinematics of the wings operating in ground effect is adjusted such 
that Pave,iGE/Pave,OGE = 100% at any time, the force ratio Fz,ave,iGE/Fz,ave,OGE is likely 
to increase. Among other factors, the feathering angle is very likely to change passively, 
when in ground effect, due to compliance of the wing [HI [271 [ig. Such effects would 
need further investigation. 

4. Conclusions 

The aerodynamic ground effect in fruitfly sized insect takeoff has been studied 
numerically using high performance computing. The three-dimensional incompressible 
Navier-Stokes equations were solved using a pseudo-spectral method with volume 
penalization using the FluSi open source code HQl HSl [H], in order to obtain the 
flow field and the aerodynamic forces acting on the insect. The takeoff trajectories 
were calculated using a simple flight dynamics solver that accounts for the body weight, 
inertia, and the legs thrust. A series of computations has been carried out to explore 
the parametric space of the model. A natural voluntary takeoff of a fruitfly, modified 
takeoffs with different kinematics and leg model parameters, and hovering flights have 
been compared. 

We found that the ground effect during the natural voluntary takeoff is negligible. 
The wingbeat averaged forces only differ by less than 1% of the weight. The aerodynamic 
power differs by less than 0.5%. 

In the modified takeoffs, we decreased the leg strength. As a consequence, the rate 
of climb decreased and the ground effect became significant. Surprisingly, the vertical 
force did not always increase. It even dropped in some of the cases that we considered. 
This is an unsteady effect related to the vortex rings bouncing off the ground surface. 
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Figure Al. Schematic drawing of the setup. 


To better understand the mechanism of the adverse ground effect, we considered 
hovering near a flat ground surface, being the limiting case of zero rate of climb. In 
that case, the fountain effect produced a large upward force on the insect’s body. The 
net ground effect was therefore positive. However, the aerodynamic force acting on the 
wings in ground effect was sometimes less than when the wings operate out of ground 
effect. The most signihcant decrease was observed during the first 15 wingbeats. Note 
that this is a much longer time period than a typical takeoff. At long time hovering, 
the effect was either positive or negative, depending on the wings kinematics. 

The parameter space in the takeoff problem is very large. In the present study, 
we focused on the legs thrust and wing kinematics. However, the aerodynamic ground 
effect may also be sensitive to the Reynolds number, because the structure of the wake 
at high Re is signihcantly different from that at low Re. Since high Reynolds number 
computations are costly, they are beyond the scope of the present study, but it is an 
important question for future research. 

Appendix A. Numerical validation of the ground plane modelling using the 
volume penalization method. 

A numerical study of the ground effect during hovering flight was carried out by Gao 
and Lu [13] in the two-dimensional approximation. In this section, we compare with 
some of their results. 

All quantities in [13] are presented in a non-dimensional form, and we follow the 
same conventions. The wing cross-section is an ellipse, as schematically shown in 
Fig. lAll Its major axis (chord length) is c = 1, and its minor axis is equal to 0.25. 
The motion of the ellipse centre in the horizontal direction (x direction) is given by 

Xcft) = Amcosiflirt/T), (A.l) 

where Am = 1.25, and we fix the dimensionless stroke time period to T = 271 Am- The 
vertical coordinate of the ellipse centre is constant in time. The angle between the major 
axis and the horizontal axis varies according to 

a{f) = ao — (Urn sin(27rt/T), 


(A.2) 
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where Oq = 90° and = 45°. 

The ground surface is horizontal {x direction), so that the distance D between the 
ellipse centre and the ground remains constant in time. In the present study, we vary 
this parameter between 1 and 6. 

The dimensionless density of the fluid is p = 1. The dimensionless kinematic 
viscosity is equal to = 10“^ yielding the Reynolds number 

Re = — = 100, where U = = i. (A.3) 

u T 

Our method solves the three-dimensional incompressible Navier-Stokes equations. 
The two-dimensional flow is therefore modelled by imposing the initial and boundary 
conditions constant in the direction perpendicular to the flow plane {z direction). The 
domain size in this direction is L 2 = 1. 

In the xy plane, the computational domain size is equal to x Ly = 12 x 12. The 
number of grid points is equal to x Ny = 512 x 512 {low resolution) or 1024 x 1024 
{high resolution). In the low resolution simulations, the grid step size Ax Ri 0.0234 is 
comparable to the lowest value reported in [13] (Ax = 0.025). 

The ground is modelled as a solid layer of width 0.2. Its top surface is at distance D 
below the centre of the ellipse. Smoothing of the penalization mask function (erf, 3Ax 
inwards and 3Ax outwards, see m) is only applied to the ellipse, not to the ground. 

The penalization parameter, both for the ellipse and for the ground, is equal 
to £ = 10“^ in the low-resolution simulations and 2.5 • 10“"^ in the high-resolution 
simulations. A ‘vorticity sponge’ forcing term m is introduced in the momentum 
equation equation in order to weaken the effect of the periodic boundary conditions in 
X. It is applied in two vertical layers of thickness equal to 32Ax, and its penalization 
parameter is equal to esponge = 0.1. 

The nearest distance to the ground in [13] is D = 1. Fig. \K2\ a'] displays the time 
evolution of the vertical force coefficient Cy, obtained by normalizing the vertical force 
Ay, 


Cy = 


A 


V 


O.SpU^cLf 

where U = 2'KAm/T. In the hgures, the time t is normalized to the stroke period T. 
After t/T = 1, oscillations of Cy are mainly described by the second harmonic. They 
become apparently periodic after t/T = 4. 

A series of simulations has been carried out in order to determine how the time 
averaged force coefficients depend on D. Their parameters correspond to the low 
resolution, as dehned above. As shown in Fig. \K2\ h']. the minimum of Cy is observed 
in our simulations at about the same D as in [13]. We conclude that, in this two- 
dimensional validation case, our results are in reasonable agreement with the reference 
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(a) (6) 




Figure A2. Results of the numerical simulations, (a) Time evolution of the 
vertical force coefficient Cy at the distance from the ground D = 1. {b) Mean vertical 
force coefficient versus distance from the ground. 



Figure Bl. Influence of the domain size in the vertical direction. Relative 
error is calculated for the aerodynamic force magnitude and the aerodynamic power, 
and plotted versus the distance from the ‘ceiling’ of the flow domain. 


Appendix B. Influence of the domain size in the vertical direction. 

In our numerical simulations of takeoffs, the insect approaches the top of the 
computational domain as time increases. The simulations have to be stopped before the 
finite domain size begins to infiuence the results significantly. To quantify this effect, 
we compare the aerodynamic forces and power in the voluntary takeoff, computed using 
two different vertical domain sizes, = 6R and = 8R. We consider the OGE case 
in which the bee is nearer to the top of the domain. 

Figure IBII shows the relative error for the aerodynamic force magnitude, 
{Fabs,ave,OGEfiR “ Fabs,ave,OGEm) / Fabs,ave,OGEm ' 100%, and for the aerodynamic power, 
{Pave,OGEfiR “ Pave,OGE, sr) / Pave,OGE,m ’ 100%. The horizontal axis shows the distance 











































Aerodynamic ground effect in fruitfly sized insect takeoff 


25 


from the insect body point of reference to the ‘ceiling’, he., the penalization layer at 
the top of the domain, normalized to the wing length R. We use the data for the 
case Lz = 8R to calculate this distance. When the domain size is equal to = QR, 

the duration of the takeoff spans four complete wingbeats, and the wings leave the 
domain during the 5th wingbeat. Hence, the 4th wingbeat is the most sensitive to the 
aerodynamic interaction with domain boundary. The maximum difference for the force 
and for the power is 0.6% and 1.6%, respectively. During the first three wingbeats, the 
error is less than 0.5%. 
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